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^ ■ ABSTRACT 

^ ! Vortex solutions in the two Higgs doublet electroweak model are constructed, and 

H ; 

5t , their stability to small perturbations is studied. The most general perturbation is 
decomposed into angular momentum modes, the least stable mode is identified, and the 
linearised energy change of the vortex under this perturbation is calculated numerically 
for various choices of parameters, thus determining whether or not the string is stable. 
It is found that, for realistic values of the Higgs mass and Weinberg angle, the string 
is unstable. 



1. Introduction 

In the last two years, interest has been generated in the study of 'embedded' defects. 
These occur when a theory which has topological defects is contained in a larger theory 
and the defects remain solutions of the field equations. General conditions for this to 
take place have been studied by Vachaspati et al . A crucial question is whether the 
defect is still stable. Generally we will not expect them to remain topologically stable 
and the stability* becomes a dynamical question. 

In this paper we consider embedded vortex solutions. A prototype example is the 
semi- local string ; the Nielsen- Olesen vortex is embedded in the semilocal model 
where the gauge group SU {2) global x U{l)iocai is spontaneously broken to U {l)giobai by a 
complex Higgs doublet. The vacuum manifold of the Higgs is now and since the first 
homotopy group of this is trivial the solution is no longer topologically stable. However 
the solutions may still be classified into topological sectors separated by infinite energy 
since finite energy requires that the Higgs field at spatial infinity lies on a C/(l) orbit on 
S^. The question is then whether the solution is dynamically stable in a given sector. 
There is one free parameter in the model, the ratio of the quartic coupling to the gauge 
coupling and the string is stable for a finite parameter range, t^'^^ 

When the SU (2) symmetry is gauged we recover the electroweak model, SU (2) l x 
U{1)y U{l)em- The vortex solution persists in the model for any value of the Wein- 
berg angle '^'^1. Now all configurations are in the vacuum sector i.e. any solution can be 
continuously deformed to the vacuum. The parameter space is now two-dimensional, 
the parameters being the ratio of the Higgs-mass to the Z-mass and the Weinberg 
angle. The string has been shown to be stable for a range in parameter space I^'^l. 
However taking the experimental lower bound on the scalar Higgs and the observed 
value of the Weinberg angle, it was shown that if the minimal standard model is the 
physically realised model then electroweak-strings are unstable. The calculation was 
refined by including finite temperature and quantum effects and while the stability 

* We are only considering the stability under small perturbations. Hence, the word 
'stability' should always be understood as 'meta-stability'. 
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of the string is improved the physical values are still outside the region of stability for 
temperatures close to the electroweak phase transition. This is somewhat disappointing 
since the prospect of a stable classical vortex is an exciting one. It could possibly be 
produced in future particle accelerators. However it is not yet understood how to cal- 
culate the cross-section for what, in quantum theory would be a coherent state of many 
particles and it may well be that the available phase space is so small as to make the 
cross-section negligible. In cosmology if the strings are meta-stable, there would be im- 
plications for physics at the electroweak phase transition. Davis and Brandenberger ^^^^ 
have suggested that they could play a role in baryogenesis. We are therefore led to 
consider extensions of the standard model where the stability region might be enlarged 
sufficiently to include the physical values. Vachaspati and Watkins ^^"^^ suggested that 
the vortex could be stabilised by the introduction of a further scalar field which would 
form a bound state with the string. 

A different possibility, the one we investigate here, is that of extending the standard 
model via the introduction of a second Higgs doublet. For simplicity we do not con- 
sider the most general such model but restrict ourselves by imposing CP-invariance 
on the Higgs potential. The model has seven parameters: the Weinberg angle, five 
Higgs couplings and the ratio of the norm of the vacuum expectation values of the two 
Higgs doublets. The stability analysis follows the same scheme as that in reference 
[8] although the introduction of the second Higgs doublet somewhat complicates the 
analysis. In section 2. we discuss the vortex solutions in the two Higgs model. In 
section 3. we consider the perturbations around the string solution and show that they 
decompose into channels labelled by angular momentum and the eigenvalues of two 
discrete symmetry operators. In section 4. we locate the least stable mode and reduce 
the analysis to solving a pair of coupled Schrodinger equations. In section 5. we discuss 
the numerical work and present our results and in 6. we give our conclusions. 
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2. Two Higgs-Doublet Electroweak Strings 

We consider the minimal extension of the Higgs sector of the standard model by 
introducing a second Higgs complex doublet. We give the most general potential subject 
only to the restrictions of gauge invariance and that it has the discrete symmetry, ($i <-> 
— $i,$2 ^ —^2) (The latter is required to ensure that flavour changing neutral 
currents are not too large.) 

$2) =iAi(|$i|' - + iA2(|$2|' - vi)^ + IXsil^if + |$2| - vf - vlf 

+ A4 (($}$i)($5$2) - |^!*2r) + As (Re($}$2) - ^^1^^2 COST/) (1) 

+ Ae ^Im($J$2) - V1V2 sin ryj 
The vacuum expectation of the Higgs fields are given by 

If sin?7 ^ then there is CP-violation in the Higgs sector. This case is considerably 
more complicated than the sin = case where, as we shall see later, the vortex solution 
is CP-invariant and the perturbations then split into channels labelled CP-even and 
CP-odd. Hence we set sin 77 = 0. Also, the Ae term is the mass term for the CP-odd 
neutral scalar, Aq, which will decouple from our analysis, so we are free to set Ae = A5. 
Then the last two terms are replaced by a single term 



Ap 



{^\^2)-ViV2 . (3) 



We are considering static, bosonic field configurations. There is no time dependence 
and we choose a gauge where the zero components of the gauge fields are set to zero. 
The energy functional is given by 



E = jd^x \G%G% + \F,^Fi^ + {D^^^)\D^^-a) + V{^r,^2) 
where a = 1,2* and 



(4) 



G«. = d,W^ - djW^ + ge'^^'^W^W^, (5) 



Barred indices always take the values 1 and 2. 



Fij = diBj - djBi, (6) 

Di^-a = di^a - ligr'^Wt^a - \ig' B^^a (7) 

where we take T" = ^ our basis for L{SU{2)) (denoting the Lie algebra of 

SU{2)) with [T«, T^] = e'^^'^T''. The Weinberg angle is given by tan Qw = g'/g- The 
semiclassical masses of the W-boson and the Z-boson are, respectively, 

= + vl), M| = \a\vl + vl), (8) 

where a — [g''^ + g'^)^ . 
There are five physical Higgs bosons, a charged pair H^^ two CP-even scalars, and 
and the CP-odd neutral scalar Aq which we have already discussed. The charged 
pair have mass rri^j^ = \4,{v1 + f|) and the two remaining Higgs scalars mix through 
the mass matrix, 

/2^;2(Ai + A3)+^;iA5 (2A3 + A5K^;2 A 

V (2A3 + A5KV2 2^;i(A2 + A3) + ^?^5; ^' 

and the physical eigenstates have masses 

^loM = \ (^11 + ^22 ± ((Mn - M22)2 + 4M22)V2) (10) 

where by convention we take ttt-Ho > '^h^- 
The time-independent field equations are 

DjGI^ = - \ig (4 T'^A^a - (A*a)V'^$«) 

AA^i =Ai(|$i|' - vl)^i + Aad^il' + |$2|' -vl- vl)^i 

(11) 

A4(|$2|^ $1 - ($5$l)$2) + A5($J$i - ViV2)^2 
AA^2 =A2(|^2|^ - vl)^2 + A3(|$l|^ + |$2|^ -vl- vl)^2 
+ A4(|$l|^ $2 - ($l$2)3'l) + A5($J$2 - ViV2)^l 

In addition, we recall the usual mixing formulae: 

Zi = cos QwW^ - sin QwBi , = sin QwW^ + cos ©VF-Bj , (12) 
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The vortex solution extremising (4) is given by : 



(13) 



$i = /i(r)e™(^;;j, $2 = /2(r)e 

where the coordinates r and 9 are polar coordinates in the xy— plane and nii and m2 
are integers. On substitution of (13) into the equations of motion they reduce to 



(14) 



(15) 



- A3 (/I - vi) h - A5 (/1/2 - ^1^2) /2 = 
/2' + 7 - (^2 - f ^ - (A2 + A3) (/I - Vl) h 

- A3 (/f - vf) h - A5 (/1/2 - ^1^2) h = 

Z"-^ + a[{rm- |Z) + (m, - |z) /|) = (16) 

where primes denote differentiation with respect to r. It is clear from the equations 
that we must have mi = 1712- This can also be seen by examining the energy functional 
and noting that to prevent a logarithmic divergence in the energy integral at infinity, 
we require that the angular dependence of the two Higgs doublets must be identical 
asymptotically. We shall restrict ourselves to the case mi = m2 = 1. The equations 
are then solved together with the boundary conditions: 

/,(0)=/,(0) = Z(0) = 0. ili^ = M-S>=^Z(^) = l (17) 

Asymptotically, for large r, writing /i = vi + Si, f2 = V2 + S2 and Z = ^ + rSs and 
linearising the equations of motion when ttihoj nT'ho < "i^^z, for example when Xijo? are 
small, we obtain the following modified Bessel equations 



/>''" + T = {:^ + Y^-i + -i))s3, (19) 

where M is the mass matrix given by (9). So, for large r, we obtain ^3 oc i^i(m2r) ~ 
r~^/^ exp{—mzr). Diagonalising the mass matrix M by writing (^^) = AiCi + A2e2, 



where 61,6*2 are the eigenvectors of M, we get decoupled equations for Ai and A2, 
yielding Ai oc Ko{mHor) and A2 oc KQ^mh^r), demonstrating that 5i and 82 are linear 
combinations of r~^/^ exp(— m^Qr) and r~^/^ exp(— m/^Qr). If or ntho is greater 
than 2mz-, the asymptotics are slightly different '■'^^1 , with ^3 source terms present in 
the Bessel equations for Ai, A2 which cannot be neglected. 
We note that for the particular value of the parameters 

ta./3.^=(^)'. (20) 

setting / = f\/v\ = /2/V2 causes equations (14) and (15) to reduce to the same 
equation: 

r + 7 - (1 - f ^ - ((Ai + Xs)v! + (A3 + X,)vl) if - l)f = 0. (21) 

In this special case, (J^) = Ai(^^^^) + A2{^''^f). We know that fi/vi = /2/W2, 
and hence 82/ 5i = tan/?, yielding A2 = 0. This shows that, despite the exponential 
dependence in A2 decaying more slowly than that in Ai, we cannot just ignore the 
more rapidly decaying Ai when studying the asymptotic behaviour of Si and 62- 

We shall examine this case in detail later since it leads to a considerable reduction in 
the complexity of the stability problem. 

3. The Harmonic Decomposition of the Perturbations 

The vortex solution given in (13) is not topologically stable. The Higgs vacuum 
is S3 X S3 and the first homotopy group of this is trivial. We are investigating the 
metastability of the vortex solution i.e. whether it is a local maximum or minimum in 
configuration space. We consider infinitesimal perturbations of the vortex configuration 
and ask if the quadratic variation in the energy is positive or negative. We write 

(/i(r)e*^'+50i)' {f2{r)e^' + S4>2) 

Z = --Z(r)e0 + SZ, A ^ a, 

r (22) 

W^T^ = W+T- + W_T+, 



where we find it convenient to use the complex notation 

W±^^ {Wi ± iW2) , T± = ^ (r^ ± ZT^) . (23) 
V 2 V 2 

The perturbations can depend on the 2— coordinate and the 2— components of the vec- 
tor fields can be non-zero also. However, since the vortex solution has translational 
invariance along the z— direction, it is easy to see that the 2;— dependence in the per- 
turbations can be ignored and the 2;— components of the gauge fields can be set to zero. 
This follows from (4) where the relevant 2;— dependent terms in the integrand are: 

iGisGia + \Fi^Fi^ + (i^3^i)^(i^3*i) + (i^3^2)^(i^3^2) (24) 

This contribution to the energy is strictly non-negative and is made to vanish by setting 
the 2— components of the gauge fields to zero and also considering the perturbations 
to be independent of the 2;— coordinate. For this reason, we shall drop all reference to 
the 2;— coordinate in the calculations below and it will be understood that the energy 
is actually the energy per unit length of the string. 

The only contribution of the photon field is the second term of the integrand in (4) 
and this gives 

i(Vxa)2 

which is positive and hence we set a = 0. 
The quadratic change in the energy may then be written as 



5E = j d^x{Xa 4>-a, SZ, W^yO 



<t>a 

SZ 



(25) 



where O, the fiuctuation operator, is a second order matrix differential operator. The 
perturbations will then decompose into channels labelled by the eigenvalues of the 
operators generating the continuous symmetries of O and by the eigenvalues of the 
discrete symmetries of O. These symmetries are the subset of the symmetries of the 
energy density which are also symmetries of the vortex configuration (13). 



It has a Zi symmetry under the action of the large XJ (2) gauge transformation given 

The perturbations of the string itself, and bZ , lie in the U = \ channel and 
and W°'T°' are in C/ = — 1. Hence these sets of perturbations decouple and we can write 

SE = SE^{6(l>a, SZ) + 6E2{Xa, W_t+). (27) 

The first term is positive since it corresponds to the stability problem of the vortex 
in the abelian Higgs case which is topologically stable as discussed in reference [15]. 

— * 

Hence we set S(/)a = 0, SZ = Q. 

The configuration is axially symmetric i.e. it is invariant under the action of the 
symmetry operator generated by the generalised angular momentum operator 

K, = L, + S, + h (28) 

where and Sg are the usual orbital and spin pieces respectively of the spatial angular 

momentum operator given explicitly by 

Q _ _ 

-C'z = , {Szh)j = -iesjkh'i- , (29) 

— * 

where 1 is the unit matrix and b is any vector field (note that Sg annihilates the Higgs 
doublets). Iz is composed of a U{1) generator, Y and an SU{2) generator, 

h = -\{Y-T^). (30) 

We therefore perform a harmonic expansion of the perturbations, 

^ (31) 

where m and n take values from — oo to +oo. It is clear that this is a decomposition 
into angular momentum modes, since 

Kzie'"^^) = me'"^^ , ^^(6*^"-^)^+) = ne^^^-^)^ . (32) 



The configuration is also CP-invariant i.e. it is invariant under the combination of re- 
flection in the x-axis and complex conjugation. Therefore the channels may be labelled 
as CP-even and CP-odd. It is straightforward to check that the U (2) gauge transforma- 
tion U = diag[z, 1] swaps the two channels and hence the two instability problems are 
exactly equivalent. Hence this supplies a reality condition and we may take X^, F„(r) 
and ^n(^) to be real. 

4. Analysis of the Least Stable Mode 

We have shown that the relevant energy variation can be decomposed into a sum over 
angular momentum modes 



(33) 



Explicitly for the m*^*^ mode we split the variation in the energy into three contribu- 
tions 

6E {Xf, Fm{r), U{r)) = Eh + E, + Ew (34) 

where E^ arises solely from the perturbation in the Higgs fields, Ew is from the per- 
turbation of the W-fields and Ec is from the interaction. 
Then 



Eh = 2n I rdr 

'0 



2^, ^,,2 



+ ^ (m+|cos2eiy^)'(|xri' + |X^|') 



+ (Ai + A3)(/i^ - vl) |xri' + (A2 + A3)(/| - vl) |X^|' 
+ A3((/i' - vl) |X-|' + (/I - vl) |xr I' + A4 |/2Xr - /iX^ 
+ 2A5(/i/2-^i^2)xrx^ 



(35) 



Ec = 27Tg / rdr 
'o 



E\Y — n I — 



- ^(/i^i + /2^2) (m + 1 - asin^ Q^Z) 



(36) 



2Frr,UlZ' + ((m - (1 - 7^)) - ^'JY + ^^^^^^ (^ + r2F^) 

(37) 



An instability is most likely to develop in the core of the string, away from the vacuum 
manifold. Considering the Higgs perturbations, if we examine Eh (35) we see that a 
Higgs condensate in the core could lead to a decrease in the potential energy. The 
Higgs perturbations can only be non- vanishing in the core for m=0, else the field is 
singular there. 

If we now consider (37) we see that the only term that can make a negative 
contribution is the first term of the integrand. The physical interpretation of this 
instability has been elucidated by Perkins ^^^^ : the magnetic field S = V x Z in the 
core of the string couples to the anomalous magnetic moment of the W-bosons and they 
acquire a negative effective mass. The potential instability is thus a W-condensate in 
the core. If we perform the Taylor expansion around the origin 

Fm = OmO + am2r'^ H , Cm = rbmo + r^bm2 H , (38) 

then we find that 

= y( "^^^ i ~ '^""f ) 1 e^(— 1)^+ + o(r^) . (39) 

^ V -IttmO Sm 6 + bmO COS 6 J 

Then the requirement that the fields be non-singular at the origin gives that they can 
only be non- vanishing there for m= 0,2. First considering the case m = 2, we then 
require F2 = ^2/^, that is they have the same sign and the contribution is therefore 
positive since Z' is positive in the core. However, for the case m = the condition is 
that Fq = —^o/r and the contribution is then negative. 

We conclude that the least stable mode is the m = mode. This confirms the 
intuitive idea that we expect the least stable mode to have = since as \Kz\ 
increases one gets an increasing centrifugal barrier. We henceforth restrict ourselves to 
the case m= 0. If we examine (35-37) we note that there are no F' terms (we drop the 
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m-script from now on). We therefore complete the squares and obtain 



5E = T^ 



dr 



f + \9\fl + fl)e 2^(1 - «sin^ ewZ)^{hX^ + /2X2) 



+ 2r2(Xi'^ + X2' ) + \{acos2ewZf{Xl + Xl) 

+ 2r2{(Ai + \Ml - vDxl + (A2 + A3)(/| - vi)x2 + A4(/2Xi - /iX2)^ 

+ A3((/l - vl)Xl + (/I - vl)Xl) + 2A5(/l/2 - ^1^2)XlX2} 

+ T(F, Xi,X2,e) 



where 



and 



and 



T(F,Xi,X2,0=^ / '-{Vp^F+ ^ 



Q = ^'(1 - 7Z) + 7Z'^ - gr^ (/iX/ - h'X, + /2X2' - /2'X2) , 



P+ = (l-7^)' + i^V(/f + /|) 



(40) 
(41) 

(42) 

(43) 



The contribution from T{F, Xi, X2, ^) is positive and we have an algebraic equation for 
F, namely 



F = - 



Q 



(44) 



so that the contribution is made to vanish. We are now left with a problem in three 
variables and we express the change in the energy in the form 



Xl 

SE[Xi,X2,^] = 271 / drr(Xi,X2,0O | X2 | , 



(45) 



where O is a 3 x 3 matrix differential operator. 

We have not completely fixed the gauge, and hence there is a gauge zero mode which 
we must extract. Perturbations of the form 



(46) 



where '0 is a real L{SU (2)) valued function and the subscript denotes the unperturbed 
fields, give an infinitesimal gauge transformation of the vortex solution (13). If we now 
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require that the gauge perturbations he in (Xi,X2)C) then we can only have 
where s(r) is any smooth function. This means that perturbations given by 

( -gfi \ 

^2 = s{r) -gf2 (48) 

are pure gauge perturbations and arc annihilated by O. Hence we need a basis for our 
physical perturbations orthogonal to this. We choose 

vi) \-v{i-^z)i \ ; 

/i(l-7Z)\ ^ ' 

+ 1 ui-iz) 



iv (fi + /I ) , 



P 

where v is a mass-scale so that the basis vectors have consistent dimensions and 



p^(i-^zr+£^{fi+fi). (50) 

Before substituting (49) into (45), to make the dependence on the parameters t = 
tan/3 = V2/V1, SIT? Qwi /3i = K/<^^ explicit we rescale the fields as follows 

R = vr, viFi = fi, V2F2 = f2, Z =\aZ , (51) 



where we now set 



1 

V = a{viV2)'^ . (52) 



Pjf. and P are now re-expressed as 

P+ = {l-2cos'^GwZf + ^R^cos^QwF^ 
P=(l-2cos^ QwZf + jCos'^QwF^ 

where 



(53) 



F'^ = -F^ + tFi. (54) 
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On substitution of (49) into (45) it was a good check on our algebra that the sq terms 
vanished. After a lot of algebra we obtain 

6E^77j dR(^s!^Ms' + s^Ns! + s^Ssj (55) 

where ' now denotes differentiation with respect to R, = (si, S2), and 

mJ. -^{i ^'^). -(t: I:)- (-) 

We give the coefficients explicitly in the appendix. Note that M is diagonal as one would 
expect from choosing orthogonal coordinates. Also, it is of interest to note that A^i2 
and 5'i2, the coefficients giving the coupling between si and S2, are both proportional 
to {F1F2 — F2FI) so that for the special choice of parameters (20), which gives Fi = -F2, 
A^i2 and S'12 are zero and the equations decouple. 

We are only assured of a complete set of eigenvectors and hence of obtaining all 
possible eigenvalues when the problem is posed in a self-adjoint form. The problem as 
it is posed in (55) is not in self-adjoint form due to the non- vanishing coefficient Ni2- 
Hence we must perform a further change of variables 

If we then choose 

^ -2M^2~ ^ ^ ' — — ' 

which we solve together with the boundary condition A(0) = then the problem (55) 
is now expressed in self-adjoint form: 

5E = TT I dR {y/'^Mu' + u^Su) , (58) 
Jo 

where trivially M and S are symmetric. Now to find the coupled Schrodinger equation 
we extremise SE subject to fixing the weighted norm of u via 

2tt / RdRu^Wu = 1, 

JQ 

1 Q 



). (59) 



where 

-A/P+ 1/P+ 
We then introduce the Lagrange multipher jj. and extremise 

/■oo 

TT / di? {nf^^^Mu! + u^^Su - 2ijRu^Wu) , (60) 
which gives the coupled Schrodinger equations 

- (My/) ' + Su = 2iiRWu. (61) 
This is equivalent to extremising 

/■oo 

TT J dR [s''^Ms' + s^Ns' + s^Ss - 2fiRs^Wsj (62) 



where 



W={1 j/pj, (63) 



which gives the coupled equations 

2(Miis;)' = + 2(5ii - 2iiR)si + 2Si2S2 

2(M224)' + {Nl2Siy = 2(^22 - 2I^R/P+)S2 + 2Si2Si. 

The boundary conditions on si and S2 are taken to be 



(64) 



sl{0) + sl{0) = 1, si(oo) = S2(oo) = 0, (65) 

which is permissible due to the linearity of the problem, and the ratio Si(0)/s2(0) will 
vary according to the eigenvalue It then follows that the stability of the string is 
dependent on since for the solutions of (64), 

poo 

5E 



poo 

= fx RdRs^Ws, (66) 
^0 



and since W and W are strictly positive definite then if fj, is negative the string is 
unstable, and if /j, is positive the string is stable. 
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5. Numerical Studies 

On making the change of variables (51), equations (14-16) become 

- tPs (F| - 1) Fi - t/35 (F1F2 - 1) F2 = 

F^' + ^-il-Z)'^-t{P, + P,){Fi-l)F, 

- j/3s (Ff - 1) F2 - jf35 (F1F2 - 1) Fi = 

Z"-^ + l{l-Z)F^^O 
which are solved together with the boundary conditions 

Fi(0) = F2(0) = Z(0) = 0, Fi(oo) = ^2(00) = Z(oo) = 1 (70) 

by relaxation. To determine the minimum eigenvalue fi, the equations (64) are then 
solved by relaxation on the same lattice. Since the couplings between si and S2 are 
small, the spectrum of eigenvalues can be split into two families; one in which s is 
mainly si, and one in which s is mainly S2- These families are studied separately by 
solving equations (64) in two separate cases, namely by imposing either si(0) = 1 or 
-§2(0) = 1, to yield two minimal eigenvalues which we will call //i and //2 respectively. 

If we set tan = 1, /?i = /?2, /?3 = /?4 = /?5 = 0, then condition (20) is satisfied and the 
si and S2 equations decouple. Writing $1 = $2 = ^/"s/S causes our energy functional 
to reduces to that of the standard electroweak model, provided we set Pi — P2 — ^P, 
where y/P is the ratio of the Higgs and Z masses. In this special case, we found the 
same region of stability for S2 perturbations as that found in references [8,9], which 
provided a useful check on the reliability of our numerics. 

It was found that fii is typically positive, giving no unstable modes made principally 
of si perturbations. Since fX2 depends only very weakly on /94, a negative fXi can easily 
be made positive by increasing the value of P4 slightly without seriously affecting //2- 
In the search for parameters which give a stable vortex solution, the main difficulty is 

1 r. 



(67) 

(68) 
(69) 



in making 112 positive. Over a wide range of parameters, it was found that the crucial 
parameter in determining 112 is 



Ho 



rrir 



Cl - C2 + ( t - ^ ) C5 



+ 4 



1/2- 



(71) 



where ci = + /Ja)/^, C2 = t(/?2 + /^a), C3 = 2/?3 + /?5, C5 = /55/2. Moving around 
the parameter space on curves of fixed mHfJ'n^z typically results in little variation in 
//2- Hence, to analyse the stability properties, we can restrict ourselves to the case 
where our energy functional reduces to that of the standard electroweak theory, thus 
reproducing the stability properties of references [8,9]. Numerical studies demonstrate 
that by tweaking the parameters in the model for fixed mass ratio (71) 112 can be 
slightly increased from the value in the electroweak model, but this increase is small. 
Some sample results are displayed in table 1. Hence, since the region of stability for 
strings in the standard electroweak model is a long way from experimentally acceptable 
values of the Weinberg angle and Higgs to Z mass ratio, the increase in the region of 
stability obtained by adding a second Higgs doublet is insufficient to lead to stable 
strings for realistic parameters. 

Table 1 



sin^ Qw 


tan/? 




/32 


/53 


/55 


n 


^J'2 


0.95 


1.0 


0.8 


0.8 


0.0 


0.0 


1.6 


-0.12* 


0.95 


1.0 


0.2 


0.2 


0.2 


0.2 


1.6 


-0.12 


0.95 


1.7 


1.17 


0.8 


-0.4 


0.5 


1.6 


-0.10 


0.95 


2.8 


0.44 


0.35 


-0.12 


0.6 


1.6 


-0.13 


0.95 


1.0 


1.6 


0.0 


0.0 


0.0 


1.6 


-0.28 


0.6 


1.0 


0.2 


0.2 


0.0 


0.0 


0.4 


-0.57* 


0.6 


1.0 


0.05 


0.05 


0.05 


0.05 


0.4 


-0.57 


0.6 


0.8 


0.04 


0.17 


0.02 


0.05 


0.4 


-0.52 


0.6 


1.0 


0.4 


0.0 


0.0 


0.0 


0.4 


-0.52 


0.6 


0.2 


0.208 


0.0 


0.0 


0.0 


0.4 


-1.86 


0.6 


2.0 


1.0 


0.0 


0.0 


0.0 


0.4 


-0.45 



In table 1, TZ = Tnjj^/rn^, and = in all cases. In the starred cases, our model 
reduces to the one doublet electroweak theory. 



1 a 



6. Conclusion 

In section 2 we demonstrated that there are vortex solutions in the CP-invariant 
two Higgs doublet electroweak theory throughout the allowable parameter range. The 
'width' of such vortices depends on the masses of the two CP-even neutral Higgs par- 
ticles and h^, and on the Z-mass. In section 3 a general perturbation of such a 
vortex was decomposed into angular momentum eigenstates, and it was seen that the 
perturbations decouple into two parts - one the same as in the two singlet abelian 
Higgs model, and a second part arising from nonzero W and Higgs perturbations or- 
thogonal to the original profiles. The CP-invariance of the model was used to decouple 
the real and imaginary parts of the perturbations. In section 4, the least stable mode 
in our expansion of the perturbations was picked out after linearisation, and a change 
of coordinates was performed to make explicit the component of this mode which was 
an infinitesimal gauge transformation. The change in energy was then expressed in 
self-adjoint form, to be analysed as a Sturm-Liouville eigenvalue problem. The numer- 
ical study of the perturbations was discussed in section 5, and it was discovered that 
the most important determining factor in the sign of the energy associated with the 
perturbation was the ratio of the masses of the and Z particles. The region of 
stability was only slightly larger than that for vortices in the one doublet electroweak 
model, and we conclude that for realistic values of the Weinberg angle and Higgs mass, 
vortex solutions in the two doublet electroweak model are unstable. 
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Appendix 

In section 4. we showed that the stabihty problem is reduced to extremising 
TT J dR [s'^Ms' + s^Ns' + s^Ss - 2iiRs^Wsj 
where ' denotes differentiation with respect to R. 



(72) 



Mil = -P+M22 = 2R 
iVi2 = 8i?(l - 2 cos^ QwZ) 



5ii =2R 



F1F2 — F[F2 



2R 

+ 



(1-2 cos^ QwR^F'^P+^) + % cos^ ©w^^ + 2i?/34F2 

R 



_.r. {F,F!,-F[F2) 
-S12 - 4i? -^-^^ 



2 cos^ QwZ' -{1-2 cos^ QwZ) — 



S22 = Ai+A2- -A's 



where 



^1 =^ cos^ {FF'f - \ cos' Qw^.FF' + \ cos' ^w^.F' 



1 F - 

+ -cos^evK-^7^ + 2cos^ 29^^(1 - 2 cos^Gvt^Z)^ 
8 RP 

- cos' Gvi^ (1 - 2 cos^ QwZ){l-2 sin' QwZ) 



RP^ 



RP^ 



2R 

+ -p2 [2 cos^ Q^Z' + (1 - 2 cos^ ewZ)P'/P] 



2n - 

+ ^{1- 2 cos' QwZy 



F2 



— COS (3 



^w^^ [{\ + i?')2cos2 QwZ' -{\-R'){l- 2cos2 QwZ)P'/Pi 

FF' r 

2cos^evK „pOD (1 + -^^) X (I + i?^)2cos2evKZ'(l -2cos2evK^) 
RP P-^ 

- (i - R')(l - 2cos2 QwZ)'P'/P 



1 Q 



2 p 

A2=^^{l-2cos''QwZf 



A3 =— ^— [cos^ QwR^FF' - 2(1 - 2 cos^ 6^/^)2 cos^ QwZ'] x 



+ 

2i?2(l - 2 cos^ eH^Z)2 + i cos^ Gvi/i^^ 



-4P- 



PP. 



+ 
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